Multistability of a two component exciton-polariton fluid 
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We study the stability of a multicomponent exciton-polariton fluid under resonant excitation 
within the linear response approximation of a generalized Gross-Pitaevskii equation. We show that, 
two spatially homogeneous and independently tunable pumping lasers produce, for the same values 
of the system parameters, up to three stable solutions. Three-stability is understood by noting 
that the cavity can be either little or highly populated and, in this second case, the largest part 
of the population lies in either one of the two components. Moreover, we discuss the different 
kinds of instabilities appearing at different pumps intensities and compare them with the case of 
one-component fluids. Finally, we show that easily tunable multistable hysteresis loops can be 
performed by the system. 
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I. INTRODUCTION 

Condensates of resonantly pumped exciton-polaritons 
in semiconductor microcavities constitute a novel and ex- 
citing system for the study of fundamental physical prop- 
erties of superfluids out of equilibrium^, and for future 
device applications 23 . Being quantum superpositions of 
light and matter they are privileged candidates for the 
realization of the next generation of optical devices 2 ^ for 
example for quantum information technologies. 

From the point of view of the fundamental physical 
properties, particularly interesting is the case of a coher- 
ently pumped polariton superfluid in presence of defects. 
Here, in contrast to the corresponding equilibrium case, 
a weak residual drag force is always present even at ex- 
tremely high polariton densities^. Nevertheless, proper- 
ties paradigmatical of an equilibrium superfluids, such as 
frictionless flow of polariton bullets 5 , quantized vortices 
and metastable persistent flow^, and the appearance and 
disappearance of Cerenkov-like waves^have been recently 
observed in coherently driven exciton-polaritons. 

In view of the potential device applications, especially 
important is the unique versatility of the polaritonic sys- 
tem, which, combined with its high non-linear properties, 
have been already demonstrated to produce parametric 
scatterin^^ and bistabilit}W. Here, the implementa- 
tion of logic operations and gates comes in a natural 
way: By manipulating the non-linear properties of the 
system using several lasers which frequencies, angles of 
incidence, and intensities can be freely varied externally. 

In this paper we investigate a new realization of a two- 
component polariton system coherently driven by two 
lasers with independently tunable frequencies, angles of 
incidence and intensities. Firstly, we study the stabil- 
ity of the two polariton components when the two laser 
intensities are varied. Differently from the case of a sin- 
gle laser pump, where the system can only be bistable, 
we disclose a rich phase diagram, where either one, two 
or three stable states can coexist at given pumping con- 
ditions. Then we suggest possible easily tunable multi- 



stable hysteresis cycles when the two pumping lasers are 
varied up and down in intensities. 

Alternatively, a multicomponent polariton fluid can be 
realized by considering the polarization degrees of free- 
dom. Multistability of different polariton spin states has 
been recently proposed theoretically 12 and confirmed ex- 
perimentally^ by the observation of three stable spin 
states for a given excitation condition. For the case of 
two polarized components, multistability in space have 
also been theoretically propose d 14 * 15 * and experimentally 
observed^. For the system with two pumping lasers pre- 
sented here, the same kind of spatial multistability is 
expected but with much more complex features. Addi- 
tionally, superimposed to the spatial multistability, in- 
terference fringes will appear due to the difference in fre- 
quency and momentum of the two pumping lasers. The 
analogy with the two-component polarized case, suggests 
that, aside from the interest in investigating multistabil- 
ity, two-component polariton condensates obtained with 
independent lasers can also be used to realize switches^ 
and memories. 

The paper is organized as follows: in section II, we 
present the model used to describe the steady state 
behavior of polaritons excited by two continuous-wave 
lasers with different frequencies, wave- vectors and in- 
tensities. The results obtained within a linear response 
framework are shown in section III. In this section we 
study the number of the possible solutions, their nature 
and possible Kerr or parametric instabilities associated 
to them. As a consequence of the presence of multiple 
stable solutions, different cycles of hysteresis can be pro- 
duced by varying, along different paths, the intensities of 
the two pumping lasers. Finally, section IV contains the 
conclusions drawn from our analysis. 



II. MODEL 

T he d ynamics of resonantly-driven microcavity polari- 
ton d 8 * 18 * can be described via a Gross-Pitaevskii equation 
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for coupled cavity (tpc) and exciton (tpx) fields general- 
ized to include decay and resonant pumping (h = 1): 
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The repulsive (gx > 0) exciton-exciton interaction in- 
duces a non-linear dynamics of the eigenmodes of the 
single polariton Hamiltonian (lower and upper polariton, 
ujlp,up0*)): 
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where G %3 = gx(\^t x \ 2 + 2|^| 2 ) with i ^ j = 1,2. 
Note that the repulsive interaction term between exci- 
tons in different states is two times larger the interaction 
term between excitons in the same mode, resulting in a 
non-uniform blue-shift. The mean-field system of equa- 
tions ([5| can have up to 9 solutions, i.e. 6 solutions more 
than in the case of one pumping laser, but, as discussed 
below, only a maximum of 3 solutions are stable. 

The dynamical stability of the two-pump-frequency 
mean-field solution can be established by adding small 
fluctuations, 



Here, we assume the cavity dispersion to be quadratic, 
u;c(k) ^c(O) + k 2 /(2m c ), with rn c = 2 x 10" 5 m (m 
is the bare electron mass), we will neglect the exciton dis- 
persion and consider the case of zero detuning at normal 
incidence, uJx(k) = co>x(0) = ujc(0). The Rabi frequency 
FIr = 5.0 [meV] and the excitonic and photonic decay 
rates, Kx = &c — 0.05 [meV] are chosen in the range of 
experimental values. 

Because of the continuous decay, a stationary state re- 
quires a continuous injection of photons. Here, we con- 
sider two continuous- wave laser fields, 
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with independently tunable frequencies uji^ and mo- 
menta ki ? 2, which can be experimentally changed by 
changing the laser angle of incidence with respect to the 
growth direction. 

We study the mean-field solutions of Eq. 
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and their stability with respect to small fluctuations 
within a linear response analysis. Substituting the ex- 
pression Q into we obtain 4 contributions, two of 
which oscillate at the main frequencies uj\ and uj2 and 
the additional two at the replica (or satellite state) fre- 
quencies uji — Auo and UJ2 + Acj, where Auj = 002 — uo\. 
Similarly to what is done in the OPO regime^ 2 ^ where 
replica states in addition to the pump signal and idler 
states are neglected, here, we consider only the terms os- 
cillating at the main frequencies uji and 002- Later, see 
Eq. (J6|, we analyse the dynamical stability of the two- 
pump-frequency solution against the weak population of 
satellite states uj{ ±uj via parametric scattering processes. 
Through the paper, we will consider only dynamically 
stable two-pump-frequency solutions. In this approxi- 
mation, we obtain the following mean-field equations for 
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where the fluctuation fields can be divided into 
particle-like and hole-like excitations 9i xc (r,t) = 



£ k [ e -«"*+* 1 "u ijri0k + e^+H 2ki - k >X x ck ]. Expanding 
Eq. up to linear terms in #i,2 x o we obtain 4 terms 
oscillating at frequencies uj\ — Auj ± uj and 002 + Auj ± cj, 
which we neglect, and 4 terms oscillating at uj^ ± uj. In 
other words, we are checking the stability of our solution, 
where only the two states with frequencies uj\^ are occu- 
pied, against the weak population of the satellite states 
Ui ± uj which can be populated by parametric scattering 
processes. The fact that we consider only linear terms 
in Ui x c k and v* x ck implies that we can obtain only 
the threshold conditions for such parametric processes, 
as well as the nature of the instability, whether of Kerr- 
type or parametric- type — see later. The equations for 
u i x ck and v* x ck can be written as an eigenvalue equa- 
tion rearranging the excitations into an 8-component vec- 
tor U T = (u lx ,ui c , v lx , v lc ,U2 x ,u 2c , v 2x ,v 2c ): 
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Here matrices L^k with i ^ j are given by 
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(5) and Ljjk are given by 
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with |^x| 2 = 2(|^ | 2 + IVCI 2 ) bein S the total exci " 
tonic density. At given values of the pumping strength 
Fi and F 2 , the solutions of the mean- field equations ([5| 
are stable if all the eight eigenvalues (LP^(k), UP^(kJ) 
of Eq. ([7]) have negative imaginary part for every value 
of the momentum k. 



III. RESULTS 

For some choices of the system parameters, we find 
that the number of stable solutions can be larger than 
one. In the case of one pumping laser, the typical S- 
like shape dependence of the polariton field intensity on 
the pump strength, also referred to as optical bistability, 
can be explained in terms of the non-linear blue-shift in- 
duced by the polariton-polariton interactio n^ 9 * 21 *. When 
the laser frequency is well above the bare lower polari- 
ton dispersion, uj p > cj£,p(kp), and the pump intensity 
increases from low values, the polariton population re- 
mains small because it is hard for the laser to inject po- 
lartions with a different energy. However, increasing the 
pump power, the blue-shift pulls the polariton energy to- 
wards resonance with the pump causing the population 
to grow superlinearly and eventually to abruptly jump to 
a high value when the pump intensity reaches a critical 
value I\. In the opposite situation, when the laser inten- 
sity is decreased from high values, the polariton energy is 
blue-detuned close to the pumping laser frequency and, 
therefore, the cavity is efficiently filled by the laser even 
at low pumping intensities. In this case the polariton 
population jumps down back at low polariton densities 
for a value I2 of the pump strength lower than I\. The 
two jumps at different values of the pump intensity cause 
therefore a hysteresis cycle. As explained below, in the 
case of two-component fluids, the situation becomes even 
richer. 

We fix both laser frequencies to be blue detuned 
with respect to the bare polariton dispersion: uji j2 = 
^Lp(ki,2) + 0.3 [meV], with ki = 0.25 /mi -1 and k2 = 
0.7 //in" 1 . We plot in panel I of Fig. [I] the phase di- 
agram showing the regions with a different number of 
stable solutions (either one, two or three) as a function 
of the two rescaled pumping intensities F{ 2 = \[QxF\,2 
[meV 3 / 2 ]. In order to understand better the structure 
of this phase diagram, we show in Fig. [2] the total ex- 
citon density gx^x^ when the pump intensities F 2 is 
kept constant at different values and F{ is varied. When 
the constant pump F 2 has a small value (Fig. [2] top left 



panel), the dependence of the population on the varying 
pump intensity F{ is similar to the one-fluid case showing 
bistability with an /S-like shape. For higher values of the 
constant pump F 2 (Fig. [2] lower left panel), the number 
of possible solutions increases but just two are found to 
be stable. Finally when F 2 is further increased the set of 
possible solutions further goes up but only a maximum 
of 3 are found to be stable. 

The coexistence of three solutions, corresponding to 
the black regions of Fig. [I] can be understood as follows: 
when the two pump intensities increase from low values, 
the polariton population is small because it's energy is 
far below the laser frequencies, resulting in one stable so- 
lution. In the opposite situation, when the intensity of 
one of the lasers decreases from high values, the polariton 
population is high and its dispersion is significantly blue 
detuned with respect to the bare one. Such blue-shift 
can be sustained by any of the two lasers, thus giving 
two additional stable solutions for the same values of the 
pump intensities. Therefore a maximum of three sta- 
ble solutions can be expected. This is also evident while 
considering the partial densities for particular values of 
the pump strength at which three stable solutions are 
present (black region of Fig. [I]), e.g. F{ = ^JgxFi = 0.05 
[meV 3 / 2 ] and F 2 = ^F 1 = 0.08 [meV 3 / 2 ]. Here, 
the solution with lower total polariton density corre- 
sponds to partial densities gx\Tpi x \ 2 — 0.009 [meV] and 
gx\^2 x \ 2 = 0-023 [meV]. The other two solutions cor- 
respond to a high value of just one of the two partial 
populations: gx\^ s \ x ^ = 0.008 [meV] and gx\^2 x \ 2 = 
0.609 [meV] in one case and gx^l^ = 0.646 [meV] 

and gx\^2 x \ 2 = O- 010 [ meV ] in tne otner - Note tnat tnis 
situation is similar to the case of two-component con- 
densates obtained with two spins. However, while in the 
spin-dependent case the two lasers pump the two spin 
populations with different intensities but at the same an- 
gle and energy, here the two pumps are independent also 
in angle and in energy. This analogy is also visible in 
Eq. ([5]). However, the difference between our system of 
equations and the spin-dependent case is that here the 
interaction between different components is twice the in- 
teraction between particles in the same component. 

In panel II of Fig. [I] we plot the phase diagram for the 
same parameters as in panel I but with the two pumping 
lasers 0.4 [meV] blue detuned with respect to the bare LP 
branch. We see that the effect of the increased detuning 
is simply to stretch the phase diagram. Since the two 
pumps are further apart from the LP branch, it is more 
difficult to inject polaritons into the cavity, and thus the 
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FIG. 2: (Color online) Stability curves of the total exciton 
density gx\^x 1 2 [meV] (red dotted curves unstable solutions, 
black lines stable solutions) for fixed pump intensities as a 
function of F{, for Fi = 0.00001 [meV 3/2] (top left), F 2 = 
0.025 [meV 3/2 ] (bottom left), and F 2 = 0.08 [meV 3/2 ] (right). 
Points Ai (Bi), i — 1,2,3 correspond to the cases when the 
lower (upper) branch of the stability curve becomes unstable 
(see Fig. [3] and Point C does not have a counterpart 

in the one-fluid case and correspond to the cases when the 
second high branch of the stability curve becomes unstable 
(see Fig. |4j). 



FIG. 1: (Color online) 2D panels: phase diagram showing the 
number of stable solutions as a function of the rescaled pump 
intensities F[ j2 — y/g^Fi^ [meV 3 / 2 ]. White, green, yellow 
and black regions correspond to respectively zero, one, two 
or three stable solutions. In panels I and II ki = 0.25 /im" 1 
and k2 = 0.7 /im -1 while in panels III and IV ki = 0.0 /im -1 
and k2 = 0.7 /im -1 . In the left panels (I and III) cji,2 = 
cJLp(ki,2) + 0.3 [meV], while in the right panels (II and IV) 
^1,2 — <^Lp(ki,2) + 0.4 [meV]. The horizontal black lines lies 
at the three fixed values of F 2 corresponding to the three pan- 
els of Figj2j while the blue diagonal line is the path used to 
plot Fig. 3D panels: plots of gx\ipx 1 2 [meV] as a func- 
tion of F[ j2 with parameters equal to panel I. Stable solutions 
with higher populations are shown in green, stable solution 
with the second higher population in yellow and third stable 
solution with lower population in black. All the solutions are 
shown in the left panel. Since the upper green branches hide 
a yellow upper branch, the right panel shows only the yellow 
and black solutions. 



need for higher pump intensities. In panel IV of Fig. [T] 
we show that a similar phase diagram can be obtained 
by changing k vector of pump 1 from 0.25 to 0.0 /im -1 . 
We observe that the multistability is quite robust with 
respect to the choice of the parameters and, therefore, it 
should be within an experimental reach. An interesting 
configuration is plotted in panel III of the same figure. 
Here a region with no stable solutions appears in the cen- 
tral part of the plot (white region). The instability of this 
region can be understood by noting that with a pump at 



k = 0.0 and just slightly blue detuned from the LP and 
a pump at k = 0.7, close to the inflection point of the 
LP, it is easy to satisfy phase matching conditions for 
parametric scattering processes. For this set of param- 
eters the system is in a configuration unstable towards 
the population of satellite states by scattering processes. 
For the other three sets of parameters, shown in the re- 
maining three panels of Fig. [I] it is also possible to find 
regions of the phase diagram where no solutions are sta- 
ble. These are the regions where the proposed solution, 
given by Eq. |4j where only the two frequency states, uj\ 
and u)2 1 are populated, is not a stable solution because 
satellite states start also to be populated — our analysis 
giving the threshold for this to happen. 

To further discuss the stability of the system with re- 
spect to small perturbations, we plot the dispersion of the 
imaginary part of the excitation eigenfrequency uj = LP^ 
for several points of the stability curves shown in Fig. [2] 
We start with the cases where the lower branch of the 
stability curve became unstable at points Ai. For very 
small values of F' 2 the imaginary part of the dispersion 
(top panel of Fig. [3]) shows two peaks for given values 
of k. One peak lies at higher value k + = 0.83 /im -1 
and one peak at lower value k_ = —0.33 /im -1 with 
k + + k_ = 2ki. This two-peaks structure is a precur- 
sor of a parametric instability due to the scattering be- 
tween two particles in the component of the condensate 
with momentum ki. This situation corresponds exactly 
to the case of one component fluids. When the pump 
intensity F' 2 is slightly increased (middle panel) we ob- 
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FIG. 3: (Color online) Dispersion of the imaginary part of 
the excitation eigenfrequency uj — LP^. The three panels 
correspond to points A\ with i = 1, 2, 3 of Fig. [2j where the 
lower part of the stability curves become unstable. In blue 
(red) the parts corresponding to the scattering of two particles 
with k = ki (k = k2). 



serve 6 other peaks appearing in the imaginary part of 
the dispersion. Two of these new peaks are such that 
k + + k_ = 2ki (blue lines) while the other four can be 
combined to identify two different scattering processes 
with k + + k_ = 2k2 (red curves). This more compli- 
cated structure of the imaginary parts of the eigenvalues 
is consistent with the fact that with two components a 
richer mechanisms of scattering might occur. When the 
intensity of F' 2 is further increased (lower panel), still 4 
different scattering may occur but, in this case, it is the 
scattering between two particles with k 2 that induce the 
instability of the system. 

For the transition from stable to unstable regions of the 
higher branch of the stability curve we plot the dispersion 
of the imaginary part of the excitation eigenfrequency 
uj = LP^ for points B { in Fig. H In analogy with the case 
of fluids with one component, Tor low intensities of pump 
2, the imaginary part of the dispersion shows a peak at 
the wavevector of pump 1, a clear precursor of a Kerr 
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FIG. 4: Dispersion of the imaginary part of the excitation 
eigenfrequency uj — LP^. The three panels correspond to 
points Bi with i = 1, 2, 3 and C of Fig. [2] where the higher 
part of the stability curve becomes unstable. 



instability. When the intensity of pump 2 is increased 
new peaks appear in the imaginary part. In the case of 
point B2 two new peaks are precursors of a parametric 
instability for the state with k = 0.25 [jim~ l ] even if the 
mechanism responsible for the instability of the solution 
is still of Kerr type. It is only when the pump intensity 
F 2 is further increased (point B3) that the two peaks at 
k- = 0.15 and fc + = 0.35 [/im _1 ] became more important 
and the mechanism of instability is of the parametric 
type. Finally in the lower right panel (corresponding to 
point C) of Fig. |4] a single peak at k = k x is the precursor 
of a Kerr-type instability that ends the region with three 
stable solutions in the right panel of Fig. [2] 

Multistability also manifests itself in a hysteresis loop 
for the populations and emission intensities obtained 
with a cycle of first increasing and later decreasing the 
pumping intensities. Here, the presence of three stable 
solutions gives more complicated loops than the ones ob- 
tained for a bistable system in a single-component po- 
lariton fluid. In order to study this aspect, we calcu- 
late the exciton emission intensity at a given frequency 
uji normalized to the total exciton emission intensity, 

a i = l^tx iVdV'ix I 2 + l^lxl 2 )' a l° n g a closed path of 
varying pumping intensities. The panels A) and B) of 
Fig. [^respectively show the hysteresis cycles of ql\ when 
the two pump intensities change along either the blue or 
the higher horizontal black line of Fig.fl] In panel A), one 
starts from a low value of F[ taking = 0.13 [meV 3 / 2 ] 
so that the population of state 2 is much higher than 
the population of state 1, i.e. Qi <C 1. Increasing F[ 
the two populations smoothly evolve until F{ ^ 0.05 
[meV 3 / 2 ] and F' 2 w 0.09 [meV 3 / 2 ]. At this point F' 2 is 
too weak to sustain high population densities in state 2 
and, therefore, the system jumps to a new stable con- 
figuration, in which the populations of both states are 
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der of 0.01 to 0.3, corresponding to a population of state 
1 being smaller but non- negligible compared to state 2. 
As F[ further increases up to 0.08 [meV 3 / 2 ], oi\ again 
jumps abruptly to values of the order of 0.95. In the re- 
verse process, F[ is decreased completing the loop. Also 
in this case, the jumps from high to low values of a\ 
are shifted to the left because of the different angles and 
energies at which the lasers are pumping the two com- 
ponent of the fluid. It is worth noting that the length 
and the height of the different plateaus of the hysteresis 
loops can be efficiently tuned by carefully choosing the 
path followed by the intensities, angles and frequencies 
of the pumping lasers. 

IV. CONCLUSION 



FIG. 5: (Color online) Hysteresis cycles of ol\ — 
l^xl 2 /(l^fxl 2 + l^| 2 ) [dimensionless] as a function of F[ for 
different values of F' 2 — stable solutions are in black, while the 
hysteresis cycle performed by the system is in brown. Panel 
A: F2 = 0.14 - F[ as in the blue line of Fig. [l] Panel B: 
F2 = 0.08 [meV 3/2 ] as in the horizontal black line of Fig.[l] 



low, i.e. ai ~ 0.5. A further increase of F{ produces a 
smooth evolution of the two populations until F[ « 0.11 
[meV 3 / 2 ] when the system jumps to a third configuration 
with a population in state 1 much higher than in state 2, 
i.e. ai « 1. When we revert the variation of the pump- 
ing intensities along the same path, the jumps to states 
corresponding to intermediate and low values of a\ are 
shifted to the left of the ones just described for increasing 
F[. The multistable hysteresis loop shown in Fig. [5|V) is 
related to the fact that the two pumping lasers are at 
different pumping angles, Iq, and pumping frequencies, 
LJi. Therefore the jumps from low to high population for 
each component appear at different values of the pump- 
ing intensities, producing the multistable behavior of ql\. 

A similar situation occurs when the system evolves 
along a path on which one of the pumping intensities 
remains constant, while the other varies (black horizon- 
tal line at F' 2 = 0.08 of panel 1 in FigfTl), as shown in 
Fig. [5j3). Starting the with F\ = 0, ol\ increases smoothly 
from zero following the lower branch up to F[ = 0.06 
[meV 3//2 ]. At this point ol\ jumps from values of the or- 



To summarize, we have studied the stability of a two- 
component exciton-polariton fluid under resonant exci- 
tation of two pumping lasers with independently tunable 
frequencies, angles of incidence and intensities. We have 
studied the effect of the detuning between the laser pump 
and the bare LP branch, and discussed the different kinds 
of instability that might occur for different values of the 
parameters. We have shown that, even though the kind 
of instabilities are the same as in the one component case 
(Kerr or parametric), here the interplay between different 
instabilities in the two components of the fluid, can pro- 
duce a much richer picture. Moreover, we have shown 
that since each component of the fluid jumps between 
stable states of its population at a different values of the 
pumping intensity, the system sustain multistable hys- 
teresis loops that can be easily modulated by changing 
the parameters of the pumping lasers. Finally, we believe 
that, due to the wide range of parameters for which the 
system is multistable, and due to the increased number 
of degrees of freedom with respect to the spin-dependent 
case, mult iat ability and hysteresis loops should be within 
an experimental reach. Therefore, the novel system of 
two-component fluid is a promising candidate for the re- 
alization of optical switches and memories. 
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